Equilibrium problem for the eigenvalues of banded block 

Toeplitz matrices 



Steven Delvaux* 



Abstract 

We consider banded block Toeplitz matrices T„ with n block rows and columns. We 
show that under certain technical assumptions, the normalized eigenvalue counting measure 
of T„ for n — oo weakly converges to one component of the unique vector of measures that 
minimizes a certain energy functional. In this way we generalize a recent result of Duits and 
Kuijlaars for the scalar case. Along the way we also obtain an equilibrium problem associated 
to an arbitrary algebraic curve, not necessarily related to a block Toeplitz matrix. 

For banded block Toeplitz matrices, there are several new phenomena that do not occur 
in the scalar case: (i) The total masses of the equilibrium measures do not necessarily form 
a simple arithmetic series but in general are obtained through a combinatorial rule; (ii) 
The limiting eigenvalue distribution may contain point masses, and there may be attracting 
point sources in the equilibrium problem; (iii) More seriously, there are examples where 
the connection between the limiting eigenvalue distribution of T„ and the solution to the 
equilibrium problem breaks down. We provide sufficient conditions guaranteeing that no 
such breakdown occurs; in particular we show this if T„ is a Hessenberg matrix. 
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1 Introduction 

Let r £ N := {1, 2, 3, . . .} and let there be given a set of r x r matrices 



fc = -a, ...,/3, 



for some a, /3 G N. These matrices are encoded by the matrix-valued Laurent polynomial (also 
called symbol) 

A(z) = A^^z-'' + . . . + Apz^ . (1.1) 
For n e N define the block Toeplitz matrix Tn{A) associated to the symbol A{z) by 

where we put Ak = if fc > /3 or fc < —a. Explicitly, 

/Ao ... \ 



TniA) 



A, 



\0 



Ar- 



A-a 



Ao J 



(1.3) 



In this paper we are interested in the limiting behavior of the eigenvalues of T„(yl) for n — > oo. 
It is known that under certain technical assumptions |21j , the eigenvalue counting measure has a 
weak limit supported on a certain curve Fq in the complex plane. An example of this phenomenon 
is shown in Figure [1] for the symbol 



Aiz) 



1 

+ z 



(1.4) 



see Bottcher-Grudsky [3] for many more illustrations of this type. 

In the case of scalar banded Toeplitz matrices, r 1, it was recently shown by Duits- 
Kuijlaars [3] that the limiting eigenvalue distribution of T„(A) satisfies a (vector) equilibrium 
problem that is constructed out of the symbol. The goal of this paper is to extend this result to 
the block case r > 1. 

Let us first review some known results in the literature, following the exposition in [6l|9]. We 
denote the eigenvalue spectrum of T„(A) by 

sp T„(A) = {A e C I det(r„(A) - A/™) = 0}, 

where in general we use 1^ to denote the identity matrix of size fc by fc. Following Schmidt- 
Spitzer [T7] , we define two limiting sets of the spectrum: we define 

liminf sp T„(A) 

to be the set of all A e C for which there exists a sequence (A„)„gN, with A„ S sp r„(A) 
converging to A. Similarly we define 



limsupsp T„(A) 
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Figure 1: Point plot in the complex plane of the eigenvalues of the banded block Tocplitz matrix 
Tn{A), n — 50, with r — 2 and symbol (|1.4p . computed in Maple with 60 digit precision. 
For — > oo the eigenvalues accumulate on a curve Fg C C which consists of six analytic arcs 
connecting the points —1, —0.42, — 1.17±0.51z, 0.28±1.58z and 1.70 (using two digits of precision). 
Three arcs are emanating with equal angles from —1 and four arcs from —0.42. The limiting 
eigenvalue distribution of Tn{A) for n — oo exists as an absolutely continuous measure /xq on Fq. 

to be the set of all A G C for which there exists a sequence (A„)„£n, with A„ G sp T„(A) having 
a subsequence converging to A. 

Under certain assumptions |21) . the above limiting sets can be described in terms of the 
algebraic equation 

= f{z, A) := det(A(z) - XU). (1.5) 

Note that each entry of the matrix A{z) — XIr is a Laurent polynomial in z, by virtue of (jl.ip . 
Hence /(z. A) = dct{A{z) — XIr) is a Laurent polynomial in z as well, and we can write it in the 
form 

/(z,A)= ^ ^(A)2^ (1.6) 

k— — q 

for certain p,q E NU {0}. The coefficients fk{X) are polynomials in A of degree at most r. More 
precisely, 

1 /• f if fc = 0, 

^'^f^-[ <r-l, ifMo! ^'-'^ 

on account of (|1.5p - (11.6p . We assume that the numbers p, q in (|1.6[) are such that the outermost 
coefficients /-<j(A) and /p(A) are not identically zero as a function of A. To avoid trivial cases 
we will always assume that 

min(p,g)>l. (1.8) 

This is justified since if mm{p,q) — 0, then det(r„(^) — A/,-„) = C'o(A)/o(A)"+" for a certain 
rational function Co (A), by Proposition 15.41 below. In that case the eigenvalues of r„(A) are 
trivially obtained. 

For any A G C with /p(A) ^ 0, we consider 

zV(z,A)= /fe(A)z^+'' (1.9) 



as a polynomial in z of degree p + g. We order its roots z = z(X) (counting multiplicities) by 
absolute value as 

< |zi(A)| < |^2(A)| < . . . < |zp+g(A)|. (1.10) 

If A is such that two or more subsequent roots in (jl.lOl) have the same absolute value, then 
we may arbitrarily label them so that (jl.lOp is satisfied. For the special values of A for which 
/p(A) = 0, the polynomial (11.91) has degree less than p + g, say p + q — j, and in that case we order 
its roots zi(A), . . . ,Zp^q-j{X) as in (jl.lOp and set Zp-|_g_j+i(A) = . . . = Zp+q(A) = oo, compare 
with We also use the latter convention if A G C is such that f{z, A) = 0. Thus in that case 
we put Zj{X) — oo ior all j = 1, . . . ,p + q. 

Each of the roots Zj(X) is finite and non-zero, except when A belongs to the set 

A := {A e C I /_,(A)/p(A) - 0}. (1.11) 

By virtue of (II. 7p . the set A has cardinality |A| < 2r — 2. In particular A is empty in the scalar 
case r = 1. 

Define the set 

To :={AeC I |z,(A)| - |z,+i(A)|}. (1.12) 

For the case of scalar banded Toeplitz matrices r = 1 it is known that Fq is a curve consisting 
of a finite number of analytic arcs and having no isolated points, and moreover the eigenvalues 
of Tn{A) accumulate on Fq in the sense that 

liminf sp Tn{A) = limsupsp r„(A) = Fq. (1.13) 

These results were shown by Schmidt and Spitzer [17]. The same authors also showed that 
the limiting eigenvalue distribution /xg of T„(A) exists as an absolutely continuous measure on 
Fq. An explicit expression for the measure fiQ was obtained by Hirschman |12] . An alternative 
expression for /xq is given by (|1.18p below with fc = 0, cf. [H]. Further results about /xq in the 
scalar case r = 1 can be found in [21 IH El [121 HH] • 

For the case of banded block Toeplitz matrices, r > 1, Widom [21 showed that the above 
results remain essentially valid, provided that the following hypotheses 112 and H3 hold true. 
The hypothesis HI is stated for further reference. 

HI. The set A in ([m]) is empty. 

H2. The set Fq in (|1.12p is a subset of C of 2-dimcnsional Lebesgue measure zero. 
H3. The set Go in (|1.14p below has finite cardinality. 
In the hypothesis H3 we define the set 

Go:={AeC\Fo|Co(A)-0}, (1.14) 

with 

Co(A) :=det f— / z^'-'' {A{z) - \Ir)-^—\ , forA£C\Fo, (1.15) 

with a in (jl.ip , and where (Tq is a counterclockwise oriented closed Jordan curve enclosing z = 
and the points Zj{\), j = I, . . . ,q, but no other roots of /(z. A) = 0. In (|1.15D the determinant 
is taken of a matrix of size ra by ra and the integral is defined entry-wise. For background. 
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generalizations and alternative representations for the function Co(A) we refer to [H [51 
(Co corresponds to the function E[ip\ in [HI HI]), see also Prop. below. 
Widom shows that under the above hypotheses H2 and H3, one has that 



liminf sp T„(A) = limsupsp Tn{A) = Tq U Gq. 



(1.16) 



It can be shown that Go is empty in the scalar case r = 1, and then (|1.16p reduces to (|1.13p . Under 
H2 and H3, Widom also observes that Hirschman's expression [12 for the limiting eigenvalue 
distribution /xo remains valid. 

If hypothesis HI fails then the limiting eigenvalue distribution of T„(A) may contain point 
masses. This is implicit in Widom [21] and will be described in detail in this paper. On the other 
hand, if H2 fails then Widom's results are not true in the stated form. Usually they remain valid 
in a modified form however, see Sections 12.31 and 16.21 below. 

The failure of hypothesis H3 is more serious, and it may cause the results to break down (see 
e.g. Section [6731) . Therefore it is important to provide sufficient conditions guaranteeing that H3 
holds true. We will provide two such conditions; in both cases H2 will hold true as well. 

Proposition 1.1. (Sufficient conditions for H2 and H3). 

(a) Suppose that the set C\ro is connected and moreover Tt^ does not have any interior points. 
Then H2 and H3 hold true. 

(h) Suppose that A{z) is the symbol of a lower Hessenberg matrix, in the sense that in the 
entry-wise expansion Tn{A) = (iij)™ZQ we have ti,j = whenever j > i + 1, i.e., all the 
entries above the first scalar superdiagonal of Tn{A) vanish. Then H2 (or more generally 
H2k below) and H3 hold true. 

Proposition 11.11 will be proved in Section [5721 Incidentally, the assumption (|1.8|) implies that 
all the entries on the first scalar superdiagonal of the Hessenberg matrix Tn{A) in Part (b) are 
non-zero. 

Finally we discuss the results of Duits-Kuijlaars [9] . These authors noticed that in addition 
to the set To in (|1.12p . an important role is played by the sets 



for /c = — (7+l,...,p— 1. In the scalar case r = 1, each set is a curve consisting of finitely 
many analytic arcs. We equip every analytic arc of F^ with an orientation and we define the 
+-side (or — side) as the side on the left (or right) of the arc when traversing it according to its 
orientation. 

Duits and Kuijlaars then define the measure 



on the curve F^. Here dX denotes the complex line element on each analytic arc of Ffe, according 
to the chosen orientation of F^. In addition, Zj-\.(\) and Zj_ (A) denote the boundary values of 
Zj(\) from the +- side and — side of F^, respectively. These boundary values exist for all but 
finitely many points. The definition (|1.18p is independent of the choice of the orientation of F^. 

In the scalar case r = 1, it is shown in [5] that the measures /i^ are the minimizers of a certain 
(vector) equilibrium problem from potential theory. Moreover, //^ is the weak limit for rt — >■ cxd 



Ffe :={AeC| |z,+fe(A)| = |z,+,+i(A)|}, 



(1.17) 




(1.18) 
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of the normalized counting measures of the kth generalized eigenvalues of the Toephtz matrix 
Tn{A). The usual eigenvalues correspond to fc = 0. 

In this paper we wish to extend these results to the block case r > 2. Instead of hypothesis 
112 we are then led to the following generalization H2fc: 

II2fc. Each set in (jl.l7p . k = —q + 1, isa subset of C of 2-dimensional Lebesgue 

measure zero. 

To the algebraic curve /(z, A) = we will associate an equilibrium problem, even when hypothe- 
ses HI and/or 112k fail. Then the equilibrium problem may contain point sources (if HI fails) or 
the definition of Tk in (|1.17p needs to be modified (if H2A: fails). 

The measure fxo will be one of the measures involved in the equilibrium problem. This 
measure will be the absolutely continuous part of the limiting eigenvalue distribution of r„ {A) , 
provided that hypothesis H3, or a suitable analogue thereof if H2 fails, holds true. In particular 
this will be the case for the two situations in Prop. 11.11 There is also an interpretation of the 
measure ^fc, A: 7^ 0, as the absolutely continuous part of the limiting distribution of the kth 
generalized eigenvalues of Tn{A), in the spirit of this will be briefly discussed in Section |31 

In the next section, we associate an equilibrium problem to an arbitrary algebraic curve 
/(z. A) = as in (|1.6p . which is not necessarily defined from a block Toeplitz matrix. In Section[3] 
we apply this to banded block Toeplitz matrices T„(A). In Section 0] we specialize our results to 
the case where T„(v4) has a scalar banded structure. Section [5] contains the proofs of our main 
results. Section [H] illustrates our results for some examples. Finally, Section [7] contains some 
concluding remarks. 

2 Equilibrium problem associated to an arbitrary algebraic 
curve 

2.1 Definitions 

In this section we show how an equilibrium problem can be associated to an arbitrary algebraic 
curve. We consider an algebraic curve which is written in the form 

/(z,A)- ^ /,(A)2'^=0, (2.1) 

k— — q 

where /a:(A), k = —q, . . . ,p are polynomials, and where p,q > 1 are such that the outermost 
polynomials f-q{X) and /p(A) are not identically zero. Note that the numbers p and q in (|2.ip 
do not have an absolute meaning; indeed by multiplying / with , j € Z, the indices p and q 
are shifted to p + j and q + j respectively. The reason why we write (|2.ip in its present form is 
because of the applications to banded block Toeplitz matrices. 
Denote 

r max deg/fc. (2.2) 

ke{~q,...,p} 

This definition of r is compatible with the one used before, by virtue of (jl.7p . 

Define the roots Zj = Zj{X), j = 1, . . . ,p + q as in (jl.lOp . and define the sets F^, fc = 
—q + 1, ... ,p — 1 as in (|1.17p . The structure of the set F^ is given by the following result. 

Lemma 2.1. (Structure of Tk). Let k G {—q+1, . . . ,p—l}. Then any point Xq G <C has an open 
neighborhood U C <C whose intersection with Tk is either empty, the singleton {Aq}, the entire 
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neighborhood U , or a finite union of analytic arcs moving from Aq to the boundary dU of the 
neighborhood U , with the arcs intersecting only at the point Aq. A similar statement holds true 
for Aq — oo provided that we consider on the Riemann sphere C C U {00} . The isolated 
points of Ffc all belong to A in (jl.lip . 

Lemma [2. II was observed for fc = by Widom [2T1 Page 312], based on the similar result for 
the scalar case by Schmidt and Spitzer [17] . The proof iov k ^ Q is exactly the same. See also 
Prop. [^?TU1 below for further information on F^. 

In addition to the set we also introduce 

f fc Ffc \ {isolated points of Ffc} (2.3) 

= cls{A e C \ A I |zg+fc(A)| = |z<,+fc+i(A)|}, 

for fc = — g + 1, . . . ,p — 1, where els denotes the closure of a subset of C. 

Our next goal is to provide an expression for the total mass of the measure /i^ in (|1.18l) . 
To this end we need some auxiliary definitions. The next definition is a variant of the so-called 
Newton polygon, see e.g. |11) . 

Definition 2.2. (The numbers rrifc). We denote by fc ruk the smallest concave function on 
{— g, . . . ,p} for which ruk > deg/fc for all k. Formally, 

mfc = max — - deg fi + - — - deg fj] , 

where the maximum is taken over all integers i,j with —q < i < k and k < j < p and with the 
equalities i — k and j = k not holding simultaneously. 

A graphical interpretation of Definition 12 .2l is as follows: consider the grid points (fc, rrik) G Z^, 
fc = — q, . . . ,p, and draw a line segment between (fc, ruk) and (fc + 1, ruk+i), for fc = —q, . . . ,p — 1. 
This then results in a curve which lies above the grid points (fc, deg/fc) £ Z^, and which is the 
'lowest' concave, piecewise linear curve with this property. 

Let us illustrate Definition 12.21 for two examples. 




Example 2.3. For the situation in Duits-Kuijlaars 191 we have deg/fc = 1 if k — and deg fk — 
if k 0. Then we easily find that 

(™fc)L-,= fo,i, l,...,-,i,o). (2.4) 

Let us illustrate this if q = Z and p = 5. In that case (deg/fc)|^_3 — (0,0,0,1,0,0,0,0,0), 
and Figure shows how to construct the concave, piecewise linear curve lying above the grid 
points (fc, deg/fc). From the figure we can then read off that (»Tifc)fc^_3 = (0, |, |, 1, |, |, |, |, 0). 
Finally, we observe that the number in (j2.4p . fc = —q+l,...,p—l, is precisely the total 
mass of the measure fj.^ in J^. 
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Example 2.4. Assume that q — 3, p — 5 and {deg fk)^^_^ = (0,1,3,3,1,0,2,0,1). Proceeding 
in a similar way as before, we find that 

(mfc)L-3= (o,|3>3,|^,2,|l^, 

as illustrated in Figure\^ 

Recall the definition of the set A in (jl.lll) , and choose an arbitrary but fixed labeling of the 
elements of this set, i.e., A =: Note that under hypothesis HI we have that L := |A| = 0, 

so in that case we can ignore all the arguments involving the numbers (A;)^]^ in what follows. 

We need the following analogue of Definition 12.21 

Definition 2.5. (The numbers tj^'^). Fix I S {1, . . . ,L}. We denote by fc i— m^''' the largest 
convex function on {— q, ■ ■ ■ ,p} for which m[,'^ < multA-Aj/fc for all k, where multA-A,/fc denotes 
the multiplicity of A — A/ as a factor of /fe(A). Formally, 

774'^= min (■^, — ^muhA-A,/i + ^ — ^muhA-A,/j ) , 

where the minimum is taken over all integers i,j with —q < i < k and k < j < p and with the 
equalities i — k and j = k not holding simultaneously. 

A graphical interpretation of Definition l2.5l is as follows: consider the grid points (fc, m^''') G Z^, 
k = —q, ■ ■ . ,p, and draw a line segment between (fc, m^''') and (fc + 1, for k = —q, . . . ,p—l. 

This then results in a curve which lies below the grid points (fc, multA-A,/fc) G and which is 
the 'highest' convex, piecewise linear curve with this property. 

Example 2.6. Assume thatq = 3, p = 5 and suppose that A; G A is such that (multA-Aj /fe)|=_3 — 
(1,1,0,2,1,0,1,1,2). Figure^then shows how to construct the convex, piecewise linear curve 
lying below the grid points (fc, multA-A,/fc)- We can then read off that 

(mi'')L-3- (1,^,0, 0,0,0, i, 1,2^. 

Definition 2.7. (The numbers fhk)- We define 

•fhk := ruk - [m)^'' + . ■ . + m^'^^ , (2.5) 

for fc = — (j'+l,...,p— 1, where and mj^''' are as in Definitions 12.21 and 12.51 
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Figure 4: Illustration for Example 12.61 



The number rhk in (|2.5p will be the total mass of the measure /i^ in p.lSp . Note that we 
defined rhk for k = — g + I, . . . ,p — 1. We could also define rhk for k — or k — p, by using the 
same definition (j2.5p . But in that case it is easy to see that m_q = fhp = 0. 

Lemma 2.8. The numbers rhk in (|2.5p satisfy 

rhk > 0, (2.6) 

for any k = —q + I, . . . ,p — 1. These inequalities are strict if p. 71) holds. 

Lemma 2.9. (Criteria for irik = 0). In Lemma \2.8[ the following statements give equivalent 
conditions for having equality in (|2.6p .• 

(i) rhk — for some k £ {—q + 1, . . . ,p — 1}, 

(a) rhk — for all k — ~q, . . . ,p, 

(Hi) f{z,X) = g{X)f{h{X)z), where g is a polynomial and h a rational function of X, and where 
f{z) — c^qZ^'^ + ■ ■ ■ + CpzP is a Laurent polynomial with coefficients Ck not depending on A. 

Lemmas l2.8l and I2.9l are proved in Section 1^751 The nonnegativity of rhk can also be deduced 
from the fact that it is the total mass of the positive measure /ifc. To avoid trivial statements, 
we will often tacitly assume that none of the equivalent conditions in Lemma 12.91 is satisfied. 

Here is an addendum to Lemma |2. II 

Proposition 2.10. (Connected components ofTk). Suppose that hypothesis H2k holds true. 
Then the number of compact, connected components ofTk is < rUk (and hence < r). Moreover, 
for each compact, connected component C ofTk, denote by fik{C) the total mass of the restriction 
of the measure fik in (jl.lSp to C, with /ifc(C) = if C is an isolated point ofTk. Then we have 
that 

HkiC) e N, if hypothesis HI holds, (2.7) 

and in general, 

AiGAnC 

where the .sum runs over all I — 1, . . . , L with A; G A fl C. 

Prop. [^TTUl will be proved in Section [575] The proposition was observed before by Widom [3TJ 
Page 315] in the special case where k — and TOq ■* ~ for all /; note that (|2.8p then reduces to 
(|2.7|) . In the scalar case r = 1 it generalizes Ullman's result [12] that Fq is connected. 
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2.2 The equilibrium problem 

Now we associate an equilibrium problem to (j2.1D . First we do this under the hypothesis H2A;. 
We will closely follow [6l [9] . For any measure /i on C define its logarithmic energy as 

= I I logl — —I dfJ.{x) dfi{y). 
J J \x~y\ 

Similarly, for any measures fi, v on C define their mutual energy as 

= / / log I ^ I dv{y). 
JJ \x~y\ 

Deflnition 2.11. We call a vector of positive measures v = {i^-q+i, ■ ■ ■ , t'p-i) admissible if I'k 
has finite logarithmic energy, Vk is supported on F^ in (|2.3p . and it has total mass Vki^k) = 'fhk 
for every /c = -g + 1, . . . ,p - 1, recall 

Definition 2.12. The energy functional J is defined by 

p-l p-2 

J{^)= X! ^i^k)- ^ I{vk,Vk+l) 
fc=-g+l k=-q+l 

+ E E (-4'li + 2m«~m»j/log-^d.,(A). (2.9) 

1=1 k=-q+l ' " 

The (vector) equilibrium problem is to minimize the energy functional (j2.9p over all admissible 
vectors of positive measures v. 

Note that the numbers — m^'lj^ + iraf' — in (12. 9[) are all nonpositive because of the 

convexity of fc i— > 

The equilibrium problem can be understood intuitively as follows, compare with [51 [5]. On 
each of the curves F^ (recall the assumption H2fc) we put charged particles with total charge 
TOfc. The particles on each curve repel each other. The particles on two consecutive curves 
attract each other, with a strength that is half as strong as the repulsion on each individual 
curve. Particles on different curves that are non-consecutive do not interact directly. Moreover, 
if 2mf^ ^ 

"i^k-i + tli^n we have an attracting external field acting on the particles on the 

curve Ffc. The external field is induced by an attracting point charge (also called sink) at A = A;. 
We refer to [15' 161 for background on equilibrium problems with external fields, and to [14] for 
vector equilibrium problems. 

Note that if hypothesis HI holds true then (12. 9|) reduces to 

k=-q+l k=-q+l 

This is the energy functional in [9] ; it also appears in the theory of Nikishin systems [14] . 
The following theorem generalizes a result in [6] and [9]. 

Theorem 2.13. (Equilibrium problem associated to an algebraic curve). Consider an algebraic 
curve as in (|2.1I) and define the sets Tk, Tk, k = —q+1, . . . ,p— 1 as in (jl.l7l) and (|2.3p . Assume 
that hypothesis H2k holds true. Then 
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(a) The vector of measures fl — defined in p.lSp is admissible, 

(h) There exist constants Ik (zS. such that 

+ E(-™fe-i + - log = (2.11) 

1=1 I 'I 

for A G r^;, k G {—(7 + 1, . . . ,p — 1}. Here we let fi-q and /ip be the zero measures. 

(c) fl — {pi'kYiZ^__q+i unique solution to the equilibrium problem in Def. \2.1^ 

Theorem 12.131 will be proved in Section [5751 Note that the equalities in Part (b) are nothing 
but the Euler-Lagrange variational conditions of the equilibrium problem, see also [5J E] . 

2.3 Roots with identically equal modulus 

Now we extend Theorem l2.13l to the case where hypothesis H2fc fails, i.e., the case where one or 
more sets Ffc have non-zero 2-dimensional Lebesgue measure in C. By Lemma 12.11 this implies 
that Ffc contains an open disk U. Inside this disk, two or more roots z[\) of f{z,X) have 
identically equal modulus as functions of A. If the disk U is disjoint from A, then we can label 
these roots so that they are analytic functions in U. The maximum modulus principle then 
implies that they are identically equal as functions of A, up to a constant factor of modulus 1. 
In this case we adapt the definition of the sets Ffc as follows: 

Ffc = {A G C \ A I |2;q+fc(A)| = \zq+k+iW\ and, possibly after relabeling the roots, the 
function Zq^k/zq+k+i takes infinitely many values on each open neighborhood U of A} 

U {A, G A I r4'^ > 0}. (2.12) 

This new definition guarantees that F^ is a curve: 
Lemma 2.14. Fix k e {—q + I, . . . ,p — 1} . 

(a) The set F^ in (j2.12p is a finite union of analytic arcs and points, with all of its isolated 
points belonging to A, and it satisfies Lemma \2.1\ 

(b) For any simply connected domain U <Z C \ (Ffc U A) , we can choose an ordering of the roots 
Zj{X) as in (jl.lOl) such that Y^jtX ^A^) o-nalytic for X U . Moreover, we can uniquely 

define the logarithmic derivative ^0^=1 1 11^=1 meromorphic function in 

C \ Ffc with poles at the points in A. 



Lemma [2. 141 is proved in Section [5741 ^ 
Due to Lemma [2.14[ we can uniquely define the measure on Tk (more precisely on F^) by 
means of (|1.18p . We have the following generalization of Theorem 12.131 

Theorem 2.15. (Equilibrium problem with roots of identically equal modulus). Consider the 
setting of Theorem \2.1S\ but assume that hypothesis H2k fails. Define the curves F^, F^ as in 
(|2.12p and (|2.3p and the measures fik as in (ll.lSp . taking into account Lemma \2.14\ Then 
Theorem \2.13\ remains valid. 

This theorem is proved in Section 15.51 
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3 The measure fiQ as the Umiting eigenvalue distribution 
of the banded block Toeplitz matrix Tn{A) 

Using the results of the previous section, we can associate a vector equihbrium problem to the 
algebraic equation f{z, A) = in (11.51) that is defined from the banded block Toeplitz matrix 
Tn{A). We want to show that the measure fiQ in the equilibrium problem is the absolutely 
continuous part of the limiting distribution of the eigenvalues of Tn{A). As discussed before, this 
will require the hypothesis H3 (or a suitable analogue thereof if H2 fails) to hold true. The next 
theorem should be compared with Widom's result [211 Theorem 6.1]. We define the normalized 
eigenvalue counting measure ^o,n of T„(A) as 

Mo,n := - h, (3.1) 

Aesp T„(A) 

where 8\ is the Dirac measure at A and each eigenvalue is counted according to its multiplicity. 

Theorem 3.1. (Limiting eigenvalue distribution of Tn{A)). Let A{z) be such that the assump- 
tions in parts (a) and/or (b) of Prov. \Ll\ are satisfied and define Fq, Gq as in (I1.12p and (jl.l4p . 
Then 

liminf sp T„(A) = limsupsp r„(A) = Tq U Go, (3.2) 

and 

Jiirn^ j <p{z) d^iQ.n{z) = j <Piz) (i^o(2) + ^m[,'V(^;) (3.3) 



1=1 

for every bounded continuous function (j) on C 

Moreover, for each A G Gq there is a positive integer j N ( more precisely, the multiplicity 
of X as a zero of Gq ) such that for every sufficiently small open disk U around A, one has 

|C/nspT„(A)| -J, (3.4) 

for all n sufficiently large, where we take into account eigenvalue multiplicities. 

Theorem 13.11 shows that the limiting eigenvalue distribution of Tn{A) for n — ^ c» consists 

of the absolutely continuous part /xq together with a point mass of mass mg'' at each A; S A, 
I — 1, . . . , L. The theorem also shows that Gq attracts isolated eigenvalues in the spectrum of 
Tn{A). The theorem will be proved in Section 

It can be checked that rrip ^ > implies A; £ Fq. The point A; can then either be an isolated 
point of Fq or it can lie on one or more analytic arcs of Fq. 

Incidentally, the occurrence of point masses at the points of {A; G A 0} can already 

be seen at the level of the finite n matrices T„(A): 

Proposition 3.2. Let A{z) in (jl.ip be the symbol of an arbitrary banded block Toeplitz matrix. 
Then there exists a constant c G M such that 

(a) For each Xi €z A, I = 1, . . . , L, we have that 

multA-A, detT„(A(z) - XIr) > m^'^n - c, for all n G N. (3.5) 

(b) We have that 

deg^detTniA{z) - XIr) <mon + c, forallneN. (3.6) 
Prop. 1321 is established in Section [ 
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The measure and the fcth generalized eigenvalues of r„(A): discussion 
Fix k G {— <? + — 1} and define the cyclic shift matrix 

S nV (3.7) 



Let A be a parameter and consider the 'shifted' symbol 

S-\A{z) - Mr) =: A^^JXjz-'^" +... + A^,(A)z^^ (3.8) 

for suitable ak, (3k S N. We may assmxie that ak, (3k are such that the coefhcients A^a^ and Ap^ 
in p.8p are not identically zero, although this will not be essential. Note that for fc = wc can 
take ao = a and (3o = (3 as in (11.11) . 

We consider the 'shifted' block Toeplitz matrix TniS^''{A{z) - Mr)). Note that for fc > 0, 
this block Toeplitz matrix is obtained from Tn{A) — Mm by skipping its first fc rows and adding 
k new rows at the bottom of the matrix, subject to the block Toeplitz structure. A similar 
description holds for fc < 0, see also 012]. 

We define the kth generalized spectrum of Tn{A) as 

spfcT„(A) = {A e C I det(T„(5-"(A(z) - A/,))) = 0}. (3.9) 

Inspired by Duits-Kuijlaars [^9, , one may hope to interpret the measure /ifc, k ^ 0, as the 
absolutely continuous part of the weak limit of the normalized counting measures of the fcth 
generalized eigenvalues of Tn{A). This limiting distribution should then also have a point mass 
of mass m^'^ at A = A;, Z = 1, . . . , L. 

It turns out that these ideas can indeed be established, provided that a suitable analogue 
II3fc of hypothesis 113 holds true. Let us define the following analogues of the objects Gq and 
Co (A) in (iri4ll - (l]l5ll : 

Gfe:={AeC\rfe|Cfe(A)=0}, (3.10) 

and 

C'fe(A) := det ( / z^-''(A(z) - Mr)-^S''—] , for A G C \ F^, (3.11) 



^ /J,i^=i,...,Qfc 

with ak in (j3.8l) . and where Uk is a counterclockwise oriented closed Jordan curve enclosing z — Q 
and the points Zj{\)^ j = 1, . . . , q + fc, but no other roots of /(z. A) — 0. In (j3.11l) the determinant 
is taken of a matrix of size ruk by rak and the integral is again defined entry-wise. 
The hypothesis H3fc now reads as follows: 

II3fc. The set Gk in (|3.10p has finite cardinality. 

Define the normalized counting measure 

1 



n 



E '5^' (3-12) 



where again each root is counted according to its multiplicity. 

Proposition 3.3. Let fc G {—q + 1, . . . ,p — 1} be such that the hypotheses H2k and H3k hold 



true. Then the statements (|3.2p - (|3.4p in Theorem \3.1\ remain true, provided that we replace 
everywhere Fq, Gq, yUo.n, MO; Wq'' and sp hy Tk, Gk, ^J'k,n, Mfe; '^i''* '^"'^ sp;. respectively. 

Unfortunately the hypothesis H3fc is very delicate to handle, and we have been unable to 
obtain sufficient conditions in the style of Prop. [TTT] for a reasonably large class of symbols A{z). 
For this reason, we will not discuss generalized eigenvalues any further in this paper. 
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4 A case study: scalar banded matrices 



In this section we specialize our results to the case where Tn{A) is a scalar banded matrix with 
non- vanishing outer diagonals. More precisely, we assume that Tn{A) is a banded block Toeplitz 
matrix as in (jl.3p . that can be written in the scalar form 



Tn{A) 



Ap) 



V 



Ap) 



\ 



rn—q—1 



(4.1) 



(k) 

where the numbers a) G C are such that 



"j ~ "i mod r' 



for alH e N U {0} and k — —q, ■ ■ ■ ,p, and with 



(4.2) 



(4.3) 



for all i = 0, . . . , r — 1. To avoid trivial cases we again assume that min(p, q) > 1. We will see in 
a moment that the notations p and q in (j4.ip are consistent with those used before in (I1.6P . 
The representations (jl.3l) and (j4.1l) are related as follows: 



(4.4) 



and the matrices A^, k — —a, . . . , in (jl.3p are obtained by taking the submatrix formed by the 
first r rows of (j4.ip and partitioning it in blocks of size r x r as follows: 



\ 



,(-■?) 



(^0 ^1 



(4.5) 



Here we add ra — q = r\q/r\ — q zero columns at the right of the matrix in the left hand side 
of (|4.5p in order to have compatible matrix dimensions. Similarly the matrices Ak, k — 0, . . . , (3 
are obtained by taking the submatrix formed by the first r columns of (|4.ip and partitioning it 
in blocks of size r x r. 

One checks that the symbol A{z) can be written as 



Aiz)= dmg{ai'\...,ai'},)S\ 

k— — q 

where S is the cyclic shift matrix in (13.71) . There is also the alternative representation 

p 



(4.6) 



DA{znD-'= ^ z'=diag(a^,'=^...,4i\)5^ 

k— — q 



(4.7) 
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where D diag(l, z, . . . ) and 



^-(/l J)- (4.8) 



One may argue that (|4.7p is more natural than (j4.6l) . in the sense that it gives the same weight 
to all the entries on the fcth scalar diagonal of the matrix (I4.ip . From this representation we 
also obtain that 

A) - det i^XIr + ^ diag(4'=\ . . . , 4i\)5^ j , (4.9) 

recaU (IT31). 



Proposition 4.1. (Structure of f). Let Tn,{A) be as in (HTT|) . T/ien f{z,\) in ca« 
6e written in the form 

f{z, A) = /-,(A)z-'' + . . . + /o(A) + . . . + /p(A)zP, (4.10) 

where all the coefficients fk{^), k = —q, . . . ,p are polynomials in A. The outermost coefficients 
take the values 

r— 1 r—1 

/_,(A) ^ = (-l)^^-!) n ai-'K /,(A) ^ /, = (-l)f('^-i) n (4-11) 

fe=o fc=o 

so hypothesis HI holds true. For general k, the degree of the polynomial fk{^) is hounded by 



, for k = -q,...,Q, 

' ^r, for k ^ 0, . . . ,p. 



Proof. Equations (|4.10p - (l4.11l) follow immediately from (14.91) . To prove (|4.12p one can use a 
combinatorial argument in the style of [HI Proof of Prop. 2.5]; a simpler proof will be obtained 
in Example 15.31 below. □ 



Corollary 4.2. Under the conditions of Prop. \4.1[ we have that 



/ w / \r> ( r 2r 2r r \ 



(4.13) 



So the total masses rhk of the measures iik form a simple arithmetic series in the same way as 
in the scalar Toeplitz case, see Example The energy functional of the equilibrium problem 
reduces to (|2.10p . 

5 Proofs 

In this section we prove our main results. 
5.1 Some preliminaries 

First wc single out some preliminaries which will be repeatedly used in the proofs. 
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Asymptotics of the roots Zj{X) 

Consider an algebraic curve f{z, A) = as in (12.11) and define the roots Zj{X) as in (jl.lOD and 
curves Ffc as in (I2.12p . Let j e {1, . . . ,p + q} be fixed. It is well-known that there exist constants 
Sj e M, e C \ {0} and kj e N such that 



z,(A)=c,A^^ (i + 0(a-i/'^^)), (5.1) 



as A — >■ cxD with A £ C\|Jj, r^, with possibly a different constant Cj for each connected component 
of C \ [Jf^ Tk in which we let A — >■ oo. Obviously, 

Sl < ... < Sp+q, (5.2) 

because of the ordering (jl.lOp of the roots zj (A) . 

The expansion (|5.ip is an instance of a Puiseux series and the next lemma is a well-known 
result for the Newton polygon. We include the proof for completeness. 

Lemma 5.1. The numbers Sj in (|5.ip are such that 

q+k 

=deg(/_g)-mfe, (5.3) 



for any k = —q, . . . ,p, with as in Definition^ 
Proof. We start from the factorization 

/p(A) 



f[zA)^'-^\{{z-z,{X)). 



By expanding this product in powers of z, we see that the coefficient fk{X) in (|2.ip is given by 

/fe(A) = (-l)^-Vp(A)5]n^^W' (5-4) 

S J£S 

where the summation runs over all subsets S C {1, . . . ,p -1- g} with |5| — p — fc, for any k G 
{—q, . . . Then we obtain 

/fe(A)//p(A)=o(^|A|^.^-^^ ) =0(|A|^^«.^^), A^oo, (5.5) 



\ s / 

for any k e {—q, ■ . ■ ,p}, where in the second step we define Sk ■— {q + k + 1, . . . ,p + q}. Hence 

p+q 

^ Sj > deg/fc - deg/p. (5.6) 

j=q+k+l 

Moreover if Sq+k < Sq+k+i then equality must hold in (j5.6p . since in that case S = Sk yields the 
unique dominant summand in the middle term of (j5.5p . In particular this holds for k = —q: 

p+q 

Sj = deg f-q - deg fp. (5.7) 
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By subtracting (15.71) from (|5.6p we then get 

q+k 

deg/-, - ^ Sj > deg/fe, (5.8) 

j=i 

with equahty if Sq+k < Sq+k+i ■ 

Denote by rhk the left hand side of (jS.Sp . Then k n- is a concave function on {— ■ ■ ■ ,p} 
by virtue of (j5.2p . From (|5.8p we see that > deg/^, with equaUty for each k for which 
2fhk > fhk-i + TOfc+i, i.e., for each k for which the concave, piecewise Unear function that 
interpolates between the grid points {k,fhk) changes slope. Then Definition 12.21 implies that 
TOfc = TTLk, which is (15. 3p . □ 

Corollary 5.2. Under the assumption (|1.7p we have that 

Zj(A)-J>0, j = l,...,(j, 
Zj(A)-J>oo, j = q + l,...,p + q, 

as A oo wit/i A G C \ IJj, F^. particular, the set Fq m (j2.12p ('or (jl.l7p ) is compact. 

Proof. The assumption (jl.7p implies that fc is a strictly increasing function on {— q, . . . , 0} 

and strictly decreasing on {0, . . . ,p}. Thus (|5.3p implies that < for j = 1, . . . , q and > 
for j = g + 1, . . . ,p + q. The result then follows from (jS.ip . □ 

Example 5.3. Lef /(z,A) he an algebraic curve as in (|4.9D a?i(i (|4.3p . T/ien 

-r/g, ^; = i,...,g, (5_^o) 

''/p, J = q + ^,---,p + q- 

Indeed, by virtue of (|5.ip and (|4.9I) w;e /ind that 

0-/(z,(A),A) - ((-l)^A^ + ciA-'?^^- +C2A^'^^)(l + o(l)), A ^ oo, 

for certain non-zero constants C\,ci. For this expression to be zero for large A we must have 
that two out of the three exponents {r, —qsj,psj} are equal and the third is smaller; this implies 
that either Sj — —r/q or Sj — r/p, for all j — 1, . . . ,p + q. The fact that Sj — ~r/q occurs 
with multiplicity q and Sj — r/p occurs with multiplicity p, is then a consequence of the relation 
j=i '^i ~ 0) recall (|5.7p and (|4.1ip . Finally, we note that (I5.10p and (j5.3D imply (j4.13p . which 
in turn leads to (I4.12p . 

Similarly to the above discussion, for any / — I, . . . , L there exist constants sj'-* G R, cj'-* e 
C \ {0} and nf^ e N such that 

z,(A) = (A ~ A,)^^" + O ((A - A0^''+^/'='') , (5.11) 

as A — > A; with A G F^, with possibly a different value of c^^ for each connected component 

of C \ IJj, Tk in which we let A — > A;. The numbers s^'' are such that 

q+k 

=multA_;,, (/_,)- (5.12) 

for any k = —q, . . . ,p and I = 1, . . . , L. 
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Widom's determinant identity 

Proposition 5.4. (Widom's determinant identity). Let A S C &e such that the solutions Zj{X) 
of the algebraic equation f{z,X) ~ in p.Sp are pairwise distinct. Then for all n sufficiently 
large we have 

detTn{A{z) - Mr) ^Y.^sW{ws{\)r+", (5.13) 
s 

where the sum is over all subsets S C {1, 2, . . . ,p + q} of cardinality \S\ = q and for each such S 
we have 

ws{X)^i-irf^,{X)l[z,{X)-' (5.14) 

jes 

and 

Cs(A) = det j z^^-'^{A{z) - XIr)-^-\ (5.15) 

where as is a counterclockwise oriented closed Jordan curve enclosing z ~Q and the points Zj{X), 
j £ S, but no other roots of f{z,X) = 0. 

Note that (|5.14p can be written alternatively as 

wsiX) = i-irf,iX)l[z,iX), 

with S :— {1,2, . . . ,p + q} \ S. This expression has maximal modulus among the subsets S of 
cardinality |5| = g if 5 = 5o with 

So = {l,...,q}. (5.16) 

For S — So the definition of Csg in (|5.15p reduces to the one of Co in (|1.15p . 

Prop. [53] was obtained in ^T, Section 6] by means of the Baxter-Schmidt formula [T]. Note 
that Prop. [53] assumes that n is sufficiently large, say n > uq, but this is no problem since [211 
Section 6, Remark 1] guarantees that the same value of no works for all A. 

Prop. assumes that the solutions of the algebraic equation f{z. A) = are pairwise distinct. 
If this assumption fails then similar determinant formulas can be obtained, by taking a suitable 
limit of (j5.13p and using continuity. This will be hinted at in Section (5.61 

For the scalar case r = 1 it is known that 

Csix) = l[z,ixy n (^.(A)-^KA))-\ 



and then Prop. \5A\ reduces to a result in [20] 



5.2 Proof of Proposition 11.11 



Proof of Provosition l 1 . lY a) . Suppose that C\ro is connected and moreover Tq does not have any 
interior points, recalling (|1.17p . From Lemma \TJ\ we immediately obtain H2. Next we establish 
H3. Eq. (15.91) implies that for |A| large enough we can take the contour (Tq in (|1.15p to be the 
unit circle. Then we easily find that 

^ [ z^--{Aiz)~XIr)-'—) =-A-i/™(l + 0(A-i)), A^oo, 

and therefore Cq(X) = (-A"^)™(1 + 0{X-^)) =^ as A ^ oo. Hypothesis H3 then follows from 
the analyticity of Co (A) in C \ Fo. □ 
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Proof of Provosition l 1 . lY b ) . Assume that T„(A) is a Hessenberg matrix. Hence by definition, 
Tn{A) has the form (14.11) with q — I and with superdiagonal entries a\~^'' ^ for all i. We will 
need some auxiliary lemmas. 

By virtue of (I4.10p (where now g = 1) we see that the equation f{z, A) = has p + 1 roots 

Zj^Zj{\), j = l,...,p+l, (5.17) 

for a certain p £ N (which is not necessarily the same p as in (14. ID ). Basic algebraic geometry 
shows that one can choose a finite union of analytic arcs F C C so that the roots ^^(A), j — 
1, . . . ,p + 1 depend analytically on A £ C, except when A S F. We will see in a moment that 
we can define F by means of (|1.17p . and the Zj(A) as in (I1.10[) : but we are not making these 
assumptions yet. 

We consider the Riemann surface TZ associated to the algebraic equation /(z, A) = 0: it is a 
branched {p + l)-sheeted covering of C, with the analytic function Zj(A) defined for A on the jth 
sheet C, j = 1, . . . ,p + 1. These functions have a cut along the appropriate arcs of F, and the 
different sheets of TZ are glued together along these arcs. 

Lemma 5.5. Let the roots Zj{X) in (j5.17l) and the Riemann surface TZ be defined as in the 
previous paragraph. Then TZ is connected. 

Proof. Take an arbitrary point (Aq, zq) G TZ. Define the set 

Z := {z E C \ there exists A G C and a continuous path in TZ from (Aq, zq) to (A, z)}. 

Then Z is a subset of C which is both open and closed. Hence it must be the entire Riemann 
sphere C. In particular it contains the value z — 0. But to z = there corresponds only A = oo 
(use (|4.10p - (|4.1ip with q — 1). Moreover, there is a unique such point (A, z) — (oo,0) on the 
Riemann surface (use (j5.9l) with q — 1). Summarizing, we see that there is a continuous path in 
TZ from (Aq, zq) to this unique reference point (oo, 0). Since this holds true for any (Aq, zq) G TZ, 
the connectedness of TZ follows. □ 

From now on we will order the roots Zj = Zj{X), j — 1, . . . ,p + 1, by increasing modulus as 
in (jl.lOD . We also define the sets F^, fc = 0, . . . ,p - 1, as in (|1.17p . 

Lemma 5.6. Each set Tk, k = 0, . . . ,p ~ 1 is a finite union of analytic arcs and points in C. 
Hence, the sets F^ can be taken as cuts for the Riemann surface TZ. 

Proof. The proof boils down to showing that F^ does not contain a (two-dimensional) open disk 
[/ C C \ A. In that case we would have two roots Zi{X) and Zj{X) that have identically equal 
modulus in U. Their ratio must be a constant of modulus one. We then obtain a contradiction by 
using the connectedness of the Riemann surface associated to f{z. A) = (Lemma 15. 5p . and the 
fact that there is only one root zi(A) that goes to zero if A goes to oo (see (|5.9p with q — 1). □ 

Since (? = 1, (|5.13p now specializes to the form 

det(T„(A) - Mm) = ^C,-i(A)(-z,(A))-"-\ (5.18) 

with 

C,-i(A) det IjA(z) - A/,.)-^^) , (5.19) 

where cr is a counterclockwise oriented closed Jordan curve enclosing z = and the point Zj{X), 
but none of the other roots Zi(A), « G {1, . . . ,p+l}, i j. Here we write Cj-i rather than Cj to be 
consistent with (|1.15l) . The function Cj_i (A) is defined for A in the domain Vj := C\ (Fj_i UFj). 
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Lemma 5.7. The function Cj_i(A) in (|5.19p . j G {1, . . . ,p + 1}, /las on/?/ isolated zeros in Vj. 

Proof. We must show that Cj_i(A) cannot be identicahy zero in any open disk in C \ lj;.rfc. 
By the fact that the Riemann surface is connected (Lemma I5.5p and analytic continuation, this 
would imply that each of the Cj_i(A), j = 1, . . . ,p+ 1, is identically zero in Vj. But then (|5.18p 
would imply that det(r„(A) — Xlm) = which is clearly a contradiction. □ 

Combining the above two lemmas, we have now established that H2fc and H3 hold true when 
Tn{A) has Hessenberg structure. This ends the proof of Proposition II . If b) . □ 

5.3 Proofs of Lemmas [23] and [2791 

Proof of Lemma \2.8l From the definition of we trivially have that 

TOfe > — ■ — deg f-q H ■ — deg fp, (5.20) 

p+q p+q 

while from the definition of to[,'^ it follows that 

ml < — ; — mu\tx-xj-q ^ ; — muhA-A,/p, (5-21) 

p+q P+q 

for any I = 1, . . . ,L. Summing (|5.2ip for all I = 1, . . . ,L and subtracting this from (|5.20p . we 
then obtain the desired inequality (j2.6p upon using that 

L L 

deg f-q = muhA-A, /-q, deg /p = ^ muhA-A, fp- (5.22) 

1=1 1=1 

Next we check the statement about the strictness of the inequality (|2.6p . From the arguments 
in the above paragraph we see that equality in (|2.6p can be achieved only if equality holds in 
both (I5.20p and (j5.2ip . For (|5.20p . this means graphically that the grid point {k,mk) lies on the 
line segment connecting (—q, deg f-q) and {p, deg fp). From the definition of the numbers 
it then follows that each of the grid points (fc,deg/fe), k = —q, . . . ,p, must lie below this line 
segment, in the sense that 

deg fk < ^ deg + ^ deg fp, (5.23) 
p+q p+q 

for all k = —q, . . . ,p. Similarly, equality in (|5.2ip implies that 

P j ^ 

muhA-A,/fe > — : — muhA-A,/-9 H : — niuhA-A,/p, (5-24) 

p+q p+q 

for all k = —q, . . . ,p and I — 1, . . . , L. 

Now if (|1.7p holds then we have deg fo — r while the right hand side of (|5.23p is at most 
r — 1. So we cannot have rhk = in that case. □ 

Proof of Lemma \2.9l First we show that (i) implies (iii). So suppose that rhk = for some 
k G {—q + 1, . . . ,p ~ 1}. As observed before, we then have the inequalities (|5.23p - (j5.24p for all 
k = —q, . . . ,p. Summing (|5.24p for all I and subtracting this from (|5.23p . we get 

L 

deg fk-J2 mult A-A, fk<0, (5.25) 
1=1 
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for any k = —q, . . . ,p, where the right hand side was simpUfied with the help of (j5.22p . On the 
other hand, we triviahy have that 

L 

deg /fe - ^ muhA-A, fk > 0. 
1=1 

So equahty holds in (j5.25p . Tracing back the argument, we must then have equality in both (I5.23P 
and (j5.24l) . Graphically this means that each of the grid points {k, degfk), k = —q, . . . ,p must 
lie on the line segment connecting {—q, deg f^q) and (p, deg/p), and similarly each of the grid 
points (fc, multA-A, /fe), k = —q, . . ■,p must lie on the line segment connecting (— g, multA-A, /-g) 
and (p, multA-Ai/p), for any I — I, . . . , L. It is easily seen that these assertions are equivalent to 
the statement in part (iii), with the rational function h{X) given by 

L 

■■= " ;^^)(muiu_,,(/p//_,))/(p+g)^ 

1=1 

So we showed that (i) implies (iii). The proof that (iii) implies (ii) can be obtained (in a simpler 
way) by reversing the above arguments. □ 

5.4 Proof of Lemma 12.141 

The proof of (a) follows again by mimicking the argument of Schmidt and Spitzer [T7]. For 
Part (b) , let [/ C C \ (Ffe U A) be a simply connected domain. For fixed X & U let ji ,j2 e N U {0} 
be such that 

\Zq+k-n-lW\ < IZq+k-j^Wl = ■■■ = \Zq+k+j2W\ < \Zq+k+j2 + lW\: 

where we set zq = and Zp+q+i = oo if necessary. Since [/ C C \ Ffc we have that either 
i2 = on U , or else ji and j2 both take a constant value on U . In the case where j2 = the 
analyticity of nj=i Zj{\) on U follows immediately from [5", Proof of Prop. 3.5]. So we can focus 
now on the case where ji and j2 are constant on U . Since J7 C C \ A, none of the roots Zj, 
i = q + k — ji,...,q + k+j2 can take the value or oo. Hence by the fact that U is simply 
connected, there exists a labeling so that each of the functions Zg+fe-^j^ (A), . . . , Zg+fc+^j (A) is 
analytic in U, with the pairwise ratios being constants of modulus 1. On the other hand, the 
argument in [51 Proof of Prop. 3.5] shows that H^ii "'^ ^ Zj{X) is analytic in U. Combining all 

these observations, we obtain the required analyticity of Hjii Zj{X) in U. 

Finally, the statement about the logarithmic derivative follows since if z(A) and z(A) are ana- 
lytic functions of X € U that are identically equal up to a constant factor, then their logarithmic 
derivatives are the same: z'(A)/z(A) = z'(A)/z(A), for all X E U. □ 



5.5 Proofs of Proposition 12.101 and Theorems 12.131 and 12.151 

In this section we prove Prop. [2TT0l and Theorems 12.131 and 12.151 The proof of Theorem 1 2 . 1 31 will 
closely follow [S] and especially [B]. 
Define the function Wk by 

q+k 

M^)^l[zj{X), AeC\rfc, (5.26) 

3 = i 
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for k — —q + 1, . . . ,p — 1. Occasionally we will also consider Wk for the indices k ~ —q oi k ~ p. 
We rewrite (|1.18p as 

d^ik{\) = TT- ^TTT 7TT d\. (5.27) 

27ri \Wk+[\) Wk-{X) J 

From Lemma 12.141 we know that w'f./wk exists as a meromorphic function on C \ with poles 
at the points of A. The following proposition gives more detailed information. 

Proposition 5.8. (The function w'f^/wk)- Let k G {— 9+ 1, . . . ,p— 1} and recall Definitions 
and \2.5\ Then the following statements hold true: 

(a) For any Aq G C \ A, there exists k G N = {1, 2, 3, . . .} such that 

Wk[\) 

as A Aq with A G C \ Ffc . We have k = \ for all hut finitely many Aq . 

(h) Near oo there exists k G N such that 

w'kW ^ deg(/_g) -ruk ^ 
Wfe(A) A 

as A ^ oo with A G C \ Ffc . 

(c) Near the point A;, I G {1, . . . , L}, there exists k G N such that 



Wk{X) A -A; 



°"'"--'^7'-^ °((A-A„-«/-). 



as A ^ A; with A G C \ F^ . 



Proof. Part (a) can be shown as in [5], for example. Now we turn to proving Part (b). Recalling 



the notation Sj in (|5.ip . we obtain 



as A — )■ oo with A G C \ F^. On account of Lemma [5.11 we then obtain part (b). The proof of 
part (c) follows in a similar way from (j5.1ip - (l5.12D . □ 

Proposition 5.9. For each k — —q+ 1, . . . ,p— 1 we have that /i^ in (j5.27p is a positive measure 
on Ffc with total mass /ifc(Ffc) — rfik- 

Proof. (Compare with Prop 3.4].) Prop . [5751 implies that the density (|5.27p is locally integrable 
around all the points in (AU{oo})nFfc, and the arguments in [9 show that fik is a positive 
measure. The statement that /ifc(Ffc) = follows from the contour deformation 

Mfe(rfe) -.^ — — — - dx 

27ri Jrfc V"'fe+(^) Wk-{X) J 
-2^/^'^^ + tl^es(^,A = A.), (5.28) 
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where C is a clockwise oriented contour surrounding Tk U A, and where Res(/i, A) denotes the 
residue of /i at A. Equation (|5.28p is vahd even if one or more points A; G A he on the curve Ffc, 
thanks to the local integrability of around these points. Applying the residue theorem once 
again, now for the exterior domain of C, we find for the first term in (j5.28p that 

1 f<W ^.J<W 



d\ = -Res[^^,X^oo\. (5.29) 
From (|5.28p - (|5.29p and the residue expressions in Prop. [5751 we then obtain 

by virtue of ([23]) and (E^H). □ 
Proposition 5.10. For each k we have that 

dfikjx) _ Wfc(A) niultA-Ai(/-g) 



\~x Wfe(A) ^ A -A, 



E . */AEC\r, (5.30) 



/ log|A - x\ dfikix) - -log|u;fe(A)| + )J(multA-A,(/-,) - m^'^loglA - Xi\ + a,., (5.31) 

1=1 

if X £ <C, for a suitable constant ak- 

Remark 5.11. As in [5], each A; G A\rfe (or A; G A) is a removable singularity for the right hand 
side of (15.301) (or (j5.3ip respectively) due to the continuity of the corresponding left hand side. 

Proof o f Provosition [5.10\ (Compare with 6, Prop 3.5].) We use the contour deformation 
dfik{x) _ 1 f 1 (w'k+i^) 



'^T^i Jf^ ^~ X \Wk+{x) Wk~{x) 
1 f 1 ,,_<W^j^l^f<ix)^^^,^ 



2'Ki Jc \ - X Wk{x) Wk{X) J^^-^i \Wk{x) 



where C is a clockwise oriented contour surrounding U AU {A}. The first term in the right hand 
side vanishes since the residue of the integrand at infinity is zero. From the residue expressions 
in Prop. [5751 we then get (|5.30p . Eq. (|5.31l) follows from this by integration, see also [5]. □ 

Finally, we can finish the proofs of Theorems 12.131 and 12.151 

Proof of Theorem[MWo')~(b). With Propositions [Hlj and jSHUj in place. Theorem |213i;a)-(b) 
now follows in the same way as in [6]. □ 
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Proof of Theorem \2.1SY c). First we show that the energy functional in (|2.9I) is bounded 

from below. To this end we rewrite J{v) as 



A:— — q+1 ^ ^ /c— — g+l ^ 

/=1 k=-q+l \ \ k / J I 



This formula is easily shown with the help of (j2.9p and (j2.5p . 

The terms in the first sum in (j5.32p are all nonnegative [Ij . For the second sum in (|5.32p , we 
observe that —irik-i + 2mk — Wfe+i — Sq+k+i — Sq+k > by virtue of (|5.3p . So these coefficients 
are all nonnegative, and they are non-zero precisely when Sg+fc+i > Sg+fc. But for such k the 
curve Tk is compact and so I{i'k) is bounded from below. Finally, for the double sum in (I5.32p 
we have that — 2m^'^ +^1!^^ = s^'^j, — s^1j._|_;^ > 0, recall (|5.11l) - (|5.12p . So these coefficients 

are all nonnegative, and they are non-zero precisely when s^^^ > s'"^}_f,_^_^. But for such k we 

have that A; ^ F^ , and then standard arguments from potential theory show that the expression 
between brackets in the double sum in (j5.32D is minimized precisely when Vk is (a constant times) 
the balayage of the Dirac point mass at A; onto the curve F^; in particular this expression is 
bounded from below as well Chapter 2]. 

Summarizing, we have now established that the energy functional J{iy) is bounded from 
below. Then the proof of Theorem I2.13f c') follows from part (b) in the same way as in [6j . □ 

Proof of Provosition \2.10\ Let us first prove (|2.8p if A n C = 0. Applying contour deformation, 
we then find that 

2mJc\Wk+W Wk~[\) J 2mJ^Wk{\) 

where 7 is the disjoint union of one or more closed Jordan curves in C \ F^. More precisely, 7 
consists of a clockwise oriented loop surrounding the outer boundary of C, and a counterclockwise 
oriented loop inside each of the 'holes' of C. Now since — (logu'fc(A))', the integral of this 

quantity over any closed Jordan curve in C \ F^ is obviously an integral multiple of 27ri. So we 
obtain (|2.8p if A n C = 0. The same argument also works if A n C 7^ provided that we take 
into account the residue from the pole of w'k/wk at each A/ e AnC, Prop. lSTST c'l. thereby noting 
that muhA-Aj-g e N U {0} C Z. 

Finally, denote by K the number of compact, connected components of F^. By summing 
(|2.8p over all such components C we get the following upper bound on K: 

c \ AiGAnc / 1=1 

by virtue of (^3)1 . □ 

Proof of Theorem \2.15\ Thanks to Lemma I2.14[ the above proof of Theorem 12.131 yields Theo- 
rem 12.151 as well. □ 



24 



5.6 Proof of Theorem 13.11 and Proposition 13.21 

Proof of Prop. \3.'A Fix A/ G A and a subset 5 C {1, 2, . . . ,p + of cardinality \S\ = q. From 
(|5.14p we have that 

ws(A) = o((A-A/)'""'*^-^'(^-')-^.'=i"^") =o((A-AO"o') , (5.33) 

as A — > A/ with A G C \ IJj, F^, where the last step follows from (|5.12p . Prop. I3.2f a') then follows 
from (|5.33p and (|5.13p . provided that there is a disk U around A/ such that for all but finitely 
many A G t/ the roots to /(z, A) = are pairwise different. But this condition is generic and the 
case where it fails follows by an easy continuity argument. 

The proof of Prop. I3.2f b) is similar. □ 

Theorem l3.1l can be obtained from Widom's determinant identity, Prop. [5^ in the same way 
as in 0!^. We outline the main steps. 

Proposition 5.12. We have that 

n-^oo A — a; . r A — X ^-^ A — A; 



uniformly on compact subsets ofC\ (Fq U Go). 

Remark 5.13. As in H], each A( G A \ (Fo U Go) is a removable singularity for the right hand side 
of (15.341) . due to the continuity of the left hand side. 

Proof of Provosition \5. 12[ As mentioned before, the dominant term in Prop. [5^ for n large is 
obtained by taking S = Sq := {1, 2, . . . , g}. Then we find in the same way as in [SJ Proof of 
Corollary 5.3] and Proof of Prop. 4.2] that 

lini / ^^"■"(^) = hm 1 V ^— = hm ^ (det T„(A(z) - A/.))' 



n—^oo 



X — X n->oo n ^ X — A, n->oo n det T„ {A(z) ~ XIr 



_ w'g^{X) _ w'pjX) ^ ^ multA-A,/-g 

-ws,iX)- u;o(A)+^ A -A; ' ^'-'"^ 

uniformly on compact subsets of C \ (Fo U Go), where the last equality in (j5.35p follows from 
(j5.26p and (|5.14p . Finally, Prop. [CTUl shows that the right hand side of (j5.35p equals the right 
hand side of ((OI)) . □ 

Proof of Theorem \3.1l From the convergence of the Cauchy transforms in Prop. 15.121 we obtain 

L 

AiO,n Mo + y^"^o'''^A, 

1=1 

in the weak-star sense, i.e., (j3.3p holds for every continuous function (f> that vanishes at infinity. 

Gerschgorin's circle theorem implies that there is a compact set K such that all the measures 
{/^o,n}n are supported in K. Therefore the assumption that vanishes at infinity is redundant 
and we obtain p.3p for all bounded continuous functions. 
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Finally, we establish the claim that Go attracts isolated eigenvalues. Let Aq G Gq and take a 
sufficiently small disk U around Aq. Then from Prop. [52] we find that 



ws„{X)-"-"dctT,,{A-XIr) = Cs„{X)+0{n, A£[/, (5.36) 

for some absolute constant c with |c| < 1. We claim that wso{X) tends to a non-zero constant 
if A — > Aq. This is obvious if Ao ^ A; if Ao = A; G A then it follows by mimicking (|5.33p and 
noting that ttiq •* = due to our assumption that Ao = A; G Go C C \ To. From (j5.36p . Hurwitz' 
theorem then implies that for all n sufficiently large, there are precisely j eigenvalues (counting 
multiplicities) of T„(j4) inside U, with j being the multiplicity of Ao as a zero of Csg = Co- 

Finally let us note that, strictly speaking, the above applications of Prop. [53] again require 
that there is a disk U around Aq such that for all but finitely many A G L/ the roots to f{z, A) = 
are pairwise different. But this constraint can again be circumvented by an easy continuity 
argument. □ 



6 Examples 

6.1 Example 1: a non-degenerate case 

We now illustrate our main results for a small-size example where each of the hypotheses HI, 
112k and H3 holds true. Consider the symbol 



Aiz) 



bi ai + ciz 

C2+a2/z 62 



(6.1) 



where we assume for convenience that each of the numbers aj, Cj, j G {1, 2}, is non-zero. Then 
the block Toeplitz matrix Tn{A) has the tridiagonal form 





/ h 

C2 


ai 
b2 


0.2 





Tn{A) = 


Cl 


bi ai 
C2 62 






I 







A little calculation shows that 



/(z,A) = -C1C22: + ((61 - A)(62 - A) - aiC2 - a2Ci) - 



aia2 



C1C2 



(z-zi(A))(z-Z2(A)), 



where as usual the roots are ordered such that |^;i(A)| < |22(A)|. We now have p — q — I and 
hence there is only one relevant set 

ro = {AGC| |zi(A)| = |Z2(A)|}. 

The coefficients Gs(A) in Prop. [Ol are labeled by index sets S C {1,2} with |5| — 1; hence 
S = {1} or 5 = {2}. It can be shown that 

1 + 

CsiX) = - 



ClC2{Zi - ZjY 
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Figure 5: Eigenvalues for the symbol (j6.ip . with the values (16.31) and n = 30, computed in 
Maple with 60 digit precision. There are 2?! ~ 60 eigenvalues in total of which 48 live on the 
vertical line segment [0:4, X3] w [—1/2 — 3. 97z, —1/2 + 3. 97i], 10 live on the horizontal line segment 
[a;i,a;2] ~ [—1.61,0.61], and the final 2 are outliers lying extremely close to A = —3 and A = 2, 
respectively, cf. 



where we put 5* = {i}, i G {1, 2}, and where j G {1, 2} is the index different from i. In particular, 
we have C5(A) — if and only if Zi{X) ~ — ai/ci. From (|6.ip . this implies in turn that 

Hence, we can only have Cs(A) = if A = 61 or A = 62- For these two special A- values, the 
second solution to f{z, A) = is Zj{X) — — a2/c2; therefore we obtain that 



Go 



{61,62}, if |ai/ci| < |a2/c2|, 
0, otherwise. 



We now turn the above discussion into a numerical example by setting 

ai = 1/2, 02 = -3, 61 =2, 62 = -3, ci =4, C2 = -3. (6.3) 

In this case, the discriminant of /(z. A) — equals 

A"* + 2X^ + 16A2 + 15A - 63/4, 

whose four roots are xi « —1.61, X2 ~ 0.61, « — 1/2 + 3. 97i and X4 ~ — 1/2 — 3.97«. These are 
the branch points of /(z. A) — 0. It turns out that the set Fq C C consists of two line segments, 
one vertically connecting the branch points 2:3 and X4 and the other one horizontally connecting 
the branch points xi and X2. The two line segments intersect at A = —1/2. 
For the values (j6.3|) . the first case in (j6.2p applies and so we have 

Go = {61,62} = {-3,2}. (6.4) 

Thus for n large, T„(A) has an isolated eigenvalue near A = —3 and near A = 2, both of 
multiplicity one. Finally, Theorem 12.131 implies that the limiting eigenvalue distribution fiQ of 
Tn{A) is precisely the equilibrium measure of the set Fq. 
These considerations are confirmed in Figure [SJ 
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Figure 6: Support Tq of the limiting eigenvalue distribution for the symbol (|6.5p . It consists of 
two closed Jordan curves in the complex A-plane and an isolated point at A = — 1. 



6.2 Example 2: a degenerate case, I 

Next we study an example where both HI and H2 fail. Consider the symbol ,21, page 321] 
Then one has that 

f{z, A) = det(A(z) - A/2) = (1 - X)z^ + z + (A^ - 1) + (-1 - X)z-\ (6.6) 

Hence hypothesis HI is violated. 

Observe that the following factorization holds, 

/(z,A)-i((l-A)z + l)(z2-A-l). (6.7) 

Hence the three roots are given by {zi(A), Z2(A), Z3(A)} = {j^, (A + 1)^^^, —(A + 1)^^^}, where 
the labeling should be taken according to increasing absolute value. Since two of the three roots 
have the same absolute value in the entire complex A-plane, hypothesis H2fc is violated as well. 

Let us first check the point sources. The set A = {A;}/!^^ (|l.llll is such that L = 2 and 
Ai = —1, A2 = 1, and the relevant data are given by 





k = -l 


A: = 


k = 1 


k = 2 


rrik 


1 


2 


3/2 


1 


(1) 


1 


1/2 








(2) 











1 


rhk 





3/2 


3/2 






From this we see that the measures /^o and /^i both have total mass tjiq = rhi — 3/2. 
Taking into account Section [^31 cf. (|2.12p . the curves Fo and Fi are defined as 

Fi = {AeC| |A + l|i/2|A-l| = l}, Fo = FiU{-l}. 
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Figure 7: Eigenvalues for n = 30 for the symbol (j6.5l) , computed in Maple with 60 digit precision. 
There are 2?! — 60 eigenvalues in total of which 15 coalesce at A = — 1. Note that the eigenvalues 
closely approximate the curve in Figure [51 



The curve Fq = Fi is plotted in Figure |6l On this curve, the measures /iq and /ii are defined 
according to the density (jl.lSp . 

The failure of H2 implies, as mentioned before, that the definitions of H3 and Gq need to 
be modified. Let us do this with an ad-hoc calculation. The coefficients Cs(A) in Prop. WM are 
labeled by index sets S C {1, 2, 3} with \S\ = 1; it can be shown that 



if 5 labels the root z(X) 



I 2(A+i)(\ta-A)(A+i)V2) . if S labels the root z(\) = ±(A + 
After some simplifications, (|5.13p then reduces to 

detT M(.) XI ) MA^-l)""^ f if n is even, 

detT„(A(.) - A/,) - (XTTka^^a^ - { if , i3 „dd. 

From this, it is easy to see that Theorem 13 . 1 1 can indeed be applied. Thus the limiting eigenvalue 
distribution of Tn{A) consists of the absolutely continuous part fig on Fq, and a point mass of 
mass 1/2 at Ai = —1. Moreover, Tn{A) does not have isolated eigenvalues for n large, neither 
for n even nor for n odd. This reproduces the result in [2T] page 321]. The comparison with the 
eigenvalues of T„ (A) for n = 30 is shown in Figure \7\ 
Finally, the energy functional (j2.9p reduces to 



lM+I{'^i)-I{i^o,i^i)~l I logp^ d^.iCA)-! log^l^ di.i(A). (6.8) 

Theorem 12.151 implies that (/io, fJ-i) is the minimizer of this functional over all pairs of measures 
(z^Oj'^i) supported on Fq = Fi, with total masses uiq — rhi ~ 3/2. The last two terms in (16. 8p 
can be interpreted as an attraction of ni towards the points A = — 1 and A = 1. 
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6.3 Example 3: a degenerate case, II 

We discuss a variant of the previous example. Consider the symbol 

The algebraic equation f{z, A) = is again given by (j6.7p . The triangularity of A{z) implies the 
following factorization for the finite n determinants: 

detT„(A(z) - A/2) = detT„(22 - 1 - A) detr„ (z^^ + 1 - A) = (-1 - A)"(l - A)". 

So the limiting eigenvalue distribution of T„(A) has a pure point spectrum with point masses at 
A = — 1 and A = 1. In particular, it is not related to the measure /zq on the set Fq in Fig. |51 
Thus Theorem 13.11 breaks down in this case. The reason is that several coefficients in Widom's 
formula are identically zero, and so H3 (actually a modification thereof since H2 fails) is not 
valid. 

It is straightforward to generalize the above idea: Whenever the symbol A{z) — XIr is block 
upper triangular, or can be reduced into block upper triangular form by means of suitable row and 
column transformations, then detT„(A(z) — A/^) factorizes into two smaller-size block Toeplitz 
determinants. For such symbols A{z), the hypotheses HI and H2fc typically hold true while H3 
and Theorem 13. II both fail. 

Finally, one may argue that the above counterexamples to Theorem 13.11 are harmless, in the 
sense that in each case the eigenvalue problem for Tn{A) can be reduced into two smaller-size 
eigenvalue problems. One may wish to construct more interesting examples for which Theo- 
rem |3IT] fails. One way to construct such examples is from the symbol 

(a{z) ai.3(z)\ 

a{z) 02,3(2) , (6.10) 
as, 1(2) 03,2(2) 03,3(2)/ 

where a(z), 0^,3(2) and 03,^(2:), i — 1,2,3, are given Laurent polynomials in z. By suitably 
fine-tuning these Laurent polynomials, and especially the exponents of their highest and lowest 
degree terms in z, one may construct symbols A{z) for which HI and H2fc hold true, H3 and 
Theorem 13.11 both fail, and for which no reduction to block upper triangular form is possible. 
We leave the details to the interested reader. 

7 Concluding remarks 

1. Generalizations. The main Theorem 13.11 was stated under the following condition: either 
C\Fo is connected and Tq does not have any interior points; or Tn{A) is a Hessenberg matrix. It 
is an open problem to generalize this theorem to other classes of banded block Toeplitz matrices. 

2. Applications. We expect that our main Theorems I2.13| / I2.15l and 13.11 may be used to 
obtain some results in the theory of multiple and matrix orthogonal polynomials on the real line. 
In fact, recently several papers appeared [S1[T31[53] which apply the results of Duits and Kuijlaars 
[5] on scalar banded Toeplitz matrices, to the context of multiple orthogonal polynomials. The 
recurrence relations of these polynomials lead to a banded Hessenberg matrix. Typically this 
matrix is not exactly Toeplitz but only asymptotically. More generally, the orthogonality weights 
may be varying with n, which leads to so-called locally Toeplitz matrices |H1[T31[53]. We anticipate 
that more applications of this type may arise in the future, possibly leading to block (rather than 
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scalar) Toeplitz matrices. A first application of this kind is given in P]. Finally, we also anticipate 
that our results could be used in the context of matrix orthogonal polynomials on the real line, 
see e.g. [3 [TO]. 
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